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ABSTRACT 

We present results of model fits to afterglow data sets of GRB 970508, GRB 980703 
and GRB 070125, characterized by long and broadband coverage. The model assumes 
synchrotron radiation (including self-absorption) from a spherical adiabatic blast wave 
and consists of analytic flux prescriptions based on numerical results. For the first time 
it combines the accuracy of hydrodynamic simulations through different stages of the 
outflow dynamics with the flexibility of simple heuristic formulas. The prescriptions 
are especially geared towards accurate description of the dynamical transition of the 
outflow from relativistic to Newtonian velocities in an arbitrary power-law density 
environment. We show that the spherical model can accurately describe the data only 
in the case of GRB 970508, for which we find a circumburst medium density n oc r~^. 
We investigate in detail the implied spectra and physical parameters of that burst. For 
the microphysics we show evidence for equipartition between the fraction of energy 
density carried by relativistic electrons and magnetic field. We also find that for the 
blast wave to be adiabatic, the fraction of electrons accelerated at the shock has to be 
smaller than 1. We present best- fit parameters for the afterglows of all three bursts, 
including uncertainties in the parameters of GRB 970508, and compare the inferred 
values to those obtained by different authors. 

Key words: hydrodynamics radiation mechanisms: non-thermal methods: statistical 
gamma-ray burst: general. 



1 INTRODUCTION 

Afterglow observations of gamma-ray bursts (GRBs) have 
provided important insight into the nature of these events. 
Some of it has been direct, for example the measurements 
of redshifts (Metzger et al. 1997), or the association of some 
bursts with supernova explosions (Hjorth et al. 2003). On 
the other hand, some has been indirect, accessible only once 
the available data are interpreted within the context of a 
physical model. The commonly used fireball model (Rees & 
Meszaros 1992; Paczynski & Rhoads 1993), for instance, is 
firmly supported by extensive modelling of afterglow obser- 
vations as synchrotron radiation originating from a decel- 
erating relativistic blast wave (Wijers et al. 1997; Waxman 
1997; Sari et al. 1998; ChevaUer & Li 2000; Panaitescu & 
Kumar 2000). 

Despite the success of the aforementioned studies in 
interpreting afterglow observations within a general frame- 
work, the values derived by independent groups for the phys- 
ical parameters of individual afterglows are often substan- 
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tially different. Such is the case for the well-studied after- 
glows of GRB 970508 and GRB 980703, for which large dif- 
ferences can be found in the derived values for blast-wave 
energy, density of the circumburst medium (CBM) and mi- 
crophysics parameters from different authors (Wijers et al. 
1997; Granot & Sari 2002; Panaitescu & Kumar 2001, 2002; 
Frail et al. 2003). The CBM density seems to be especially 
unconstrained, as differences of many orders of magnitude 
can be found in the literature. One of the most important 
parameters of GRB outflows, that directly affects the in- 
ferred energetics and rate of these events, is the opening 
angle of the jet. Specifically, jetted instead of spherical out- 
flows would significantly aleviate the energy requirements 
and boost the event rate of GRBs. The first strong infer- 
ence of their presence (Harrison et al. 1999) was perceived 
as evidence for the ubiquitous role they play in the GRB phe- 
nomenon. Accumulating observations, however, have failed 
to fully confirm this picture, with many afterglows not show- 
ing any steepening in the light curves that can be attributed 
to a jet break (e.g. Racusin et al. 2009), rendering the in- 
fluence of collimation on GRB outflows for the most part 
ambiguous. All these uncertainties on the inferred physical 
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parameters of GRB blast-waves have called for refinement 
and greater precision in the methods that underlie afterglow 
modelling. 

Theoretical afterglow calculations have been continu- 
ously improved to include more precise methods of calculat- 
ing the dynamics and spectra of the source (e.g. Kobayashi 
et al. 1999; Huang et al. 1999; Rhoads 1999; Granot & Sari 
2002; Granot & Piran 2012; Pe'cr 2012). Many recent stud- 
ies (e.g. Meliani et al. 2007; Zhang & MacFadyen 2009; van 
Eerten et al. 2010a; Wygoda et al. 2011; De CoUe et al. 
2012a; van Eerten et al. 2012) arc based on high-resolution 
relativistic hydrodynamic (RHD) simulations which are es- 
sential to understand critical aspects of the outflow's dy- 
namics, like lateral spreading of jets and the transition to 
the non-relativistic phase. This allows in principle for accu- 
rate determination of spectra and light curves from simula- 
tion runs. However, this method is not suitable for iterative 
fitting of model parameters to observations due to the lim- 
itations posed by the necessary performance of numerous 
time-consuming RHD simulations. 

Recently (van Eerten et al. 2012; Leventis et al. 2012; 
see also van Eerten & MacFadyen 2012b) a new approach 
has been developed for the calculation of spectra and light 
curves that retains the accuracy of the numerical techniques, 
without requiring the long run times of simulations. While 
the methods of these studies differ, they are common in how 
they are based on sets of blast-wave simulations that span 
the parameter space. In the case of van Eerten et al. (2012) 
dynamical results of 2D simulations have been tabulated 
allowing the user to perform a straightforward immerical 
calculation of the afterglow radiation for any combination 
of the physical parameters within the explored range. Even 
so, this calculation can be lengthy and is best executed on 
a parallel computer network. 

The method of Leventis et al. (2012) is based on ID 
RHD simulations that span the entire range of dynamics, 
from ultrarelativistic to Newtonian velocities. These sim- 
ulations, however, do not account for jet features as they 
rely on the assumption of spherical symmetry. Several runs 
have been used to calibrate analytically derived scalings of 
observed synchrotron spectra. The resulting formulas have 
the unique advantage of combining the accuracy of high- 
resolution trans-relativistic simulations with the versatility 
of analytic equations. The fact that they cover a secjuenco 
of dynamical phases has motivated us to use them in order 
to fit model parameters to observational data for afterglows 
with extensive monitoring. The bursts we are mainly con- 
cerned with in this paper are GRB 970508, GRB 980703 and 
GRB 070125, all monitored in several bands from radio to 
X-ray frequencies and covering observer times from hours to 
several months. The two former are among the most studied 
afterglows with several groups publishing results they have 
obtained through afterglow modelling. 

In this work we present fit results for the afterglows of 
these bursts and investigate the extent to which a spheri- 
cal outflow can provide an adequate description of the data. 
These results also serve as a basis for comparison to model 
fits based on 2D simulations. Furthermore, the prescriptions 
of Leventis et al. (2012) enable us to examine the density 
structure of the burster's immediate environment, as a con- 
tinuous range of values for the slope of the CBM density is 
allowed. The resulting slope can then reveal unusual density 



distributions of the CBM, or confirm previous claims based 
on models with only preset structures available, typically 
constant density or a profile corresponding to a stellar-wind 
environment (oc r~^). 

The paper is organized as follows. A description of the 
observational data that have been used during fitting is pre- 
sented in Section 2. In Section 3 we illustrate the main fea- 
tures of the physical model we have used and in Section 4 
we present our main results. In Section 5 we focus on the 
inferred parameters of GRB 970508 for which we obtain the 
most reliable results. In Section 6 we discuss the implications 
of this work on afterglow physics and modelling. Finally, in 
Section 7 we conclude by summarizing our main findings. 



2 DATA 

In this study we focus on three sources: GRB 970508, 
GRB 980703 and GRB 070125. All three have well-sampled 
afterglows across the electromagnetic spectrum. In partic- 
ular they are among the few GRBs that have detections 
in multiple radio bands at hundreds of days after the initial 
gamma-ray trigger. This allows us to model the full evolution 
of the GRB blast wave from the ultrarelativistic to the non- 
relativistic phase. Another burst with afterglow monitoring 
spanning almost a decade in the radio is GRB 030329. We 
have not fit that data set as it is clear from the light curves 
that a jetted model is needed to interpret the observations 
(see van der Horst et al. 2008 and references therein). 

Since the lauch of the Swift satellite, it has become clear 
that the early (10"^ — 10^ s) afterglow behaviour of many 
bursts cannot be explained by standard afterglow models 
(Nousek et al. 2006). Energy injection into the blast wave 
has been proposed to explain the typically shallow decay 
that the optical and X-ray light curves show (e.g. Gra- 
not & Kumar 2006; Nousek et al. 2006; Zhang et al. 2006; 
Panaitescu & Vestrand 2011). Other plausible explanations 
are evolution of the shock microphysics parameters (Gra- 
not et al. 2006), or viewing angle effects (Eicliler & Gra- 
not 2006). In our sample, GRB 970508 and GRB 070125 dis- 
play an atypical behaviour, lasting in both cases up to 1.5 
days. Especially for GRB 970508, the fast-rising optical light 
curves before 1.5 days may reveal a refreshed shock, occur- 
ing when a slow shell catches up with the afterglow shock at 
later times (Kumar & Piran 2000; Granot et al. 2003). After 
1.5 days the light curves are compatible with the canonical 
afterglow decay. Processes like energy injection, refreshed 
shocks and effects due to off-axis viewing angle cannot be 
accounted for in the model we are using in this work. For 
this reason we have excluded data before 1.5 days from the 
fitted data sets of both GRB 970508 and GRB 070125. 

For GRB 970508 radio observations were performed at 
1.43, 4.86 and 8.46 GHz (Galama et al. 1998b; Frail et al. 
2000) . Near-infrared and optical data have been published at 
6 observing bands (Chary et al. 1998; Galama et al. 1998a; 
Sokolov et al. 1998; Sahu et al. 1997; Garcia et al. 1998). 
The magnitudes of the underlying host galaxy in the B, V, 
Rc and Jc bands have been presented in Zharikov & Sokolov 
(1999), while the observations in the K and U bands are suf- 
ficiently early that they are not affected by the host galaxy 
brightness. We have corrected the observed optical magni- 
tudes for galactic extinction, subtracted the host galaxy flux. 
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and converted them to fluxes. The afterglow was observed 
in X rays with BeppoSAX (Piro et al. 1998), for which we 
have converted the X-ray count rates to fluxes by assuming 
a spectral index of —1.1 over the observing band. 

GRB 980703 was observed at the same radio frequencies 
as GRB 970508 (Berger et al. 2001; Frail et al. 2003). We 
have used all the near-infrared and optical data in the H, J, 
I, R, V and B bands (Bloom et al. 1998; Castro-Tirado et al. 
1999; Vreeswijk et al. 1999). We have corrected the observed 
magnitudes for galactic extinction, but also for extinction in 
the host galaxy with E{B -V) = 0.29 (Starling et al. 2007; 
Starling 2008). The host galaxy of GRB 980703 was bright, 
not only in the optical (Frail et al. 2003) but also at radio 
wavelengths (Berger et al. 2001), and we have subtracted 
the host galaxy flux at all these wavelengths from our mea- 
sured fluxes. The afterglow has also been detected at X-ray 
energies (Vreeswijk et al. 1999), for which we have used the 
same conversion method as in the case of GRB 970508. 

For GRB 070125 we have used all the broadband data 
presented in De Cia et al. (2011). Radio observations were 
performed at 4.86, 8.46, 15 and 22.5 GHz, while millime- 
tre observations were carried out at 95 and 250 GHz. The 
data set is supplemented by observations at 12 more bands 
ranging from the near infrared to X-ray energies, including 
optical and ultraviolet bands. 

To carry out the modelling we have adopted the fol- 
lowing cosmology: Om = 0.27, Oa = 0.73 and the Hubble- 
parameter Ho = 71 kms"^ Mpc"^ so for the GRB 970508 
redshift of z = 0.835 (Metzger et al. 1997) the luminos- 
ity distance is dL = 1.64 • 10^* cm, for GRB 980703 the 
redshift z — 0.966 (Djorgovski et al. 1998) corresponds 
to dL = 1.96 • lO^^cm, and for GRB 070125 the redshift 
z = 1.547 (Cenko et al. 2008) implies that dt = 3.53-10^* cm. 

During the flt process, no data were excluded based on 
flux values. However, in the figures we present, data that are 
not significant at the 2a level are depicted as upper limits, 
for display purposes. 

2.1 Scatter of the data 

A noticeable feature of radio data is the high degree of scat- 
ter they show, especially compared to the size of the error 
bars (see Section 4). In the case of GRB 970508, notable scat- 
ter is also present in near-infrared and optical frequencies. 
In these bands it is presumably caused by the use of data 
from various telescopes without the performance of cross- 
calibration analysis. 

In the radio, interstellar scintillation affects the flux lev- 
els (Goodman 1997). Its strength diminishes as the angular 
size of the source grows and this has been used to infer the 
radius of GRB outflows (Frail et al. 1997; Taylor et al. 1997; 
Waxman et al. 1998; FraU et al. 2000; Yost et al. 2003). 
Various other groups have accounted for the effect of scin- 
tillation (Panaitescu & Kumar 2002; Chandra et al. 2008), 
especially in early data, by effectively increasing the size of 
the error bars, more so in early observer times. In our study 
we have not included the effect of scintillation in the model. 
One reason is that it does not affect the best-fit values of 
our results significantly, as the central values of the mea- 
surements are not perturbed. Another reason is the lack of 
detailed measurements for the amount of scattering material 
off the galactic plane, which makes the effect of scintillation 



in models of extragalactic sources uncertain (Chandra et al. 
2008). 



3 THE MODEL 

3.1 General description 

The model we have used is a direct implementation of the 
method presented in Leventis et al. (2012). In that pa- 
per we present simulation-calibrated flux prescriptions of 
synchrotron radiation, including self-absorption, through- 
out the entire dynamical evolution of GRB afterglows. The 
model assumes an initially ultrarelativistic spherical blsist 
wave expanding adiabatically inside a medium with a den- 
sity profile described by a power law: n{r) ocr~*^. The en- 
ergy distribution of the electrons accelerated at the forward 
shock is also assumed to be a power law. The minimum 
Lorentz factor of that distribution is calculated through the 
energy density and mass density of the shocked gas. The syn- 
chrotron spectrum is then determined through the emissiv- 
ity and absorption coefficient of these relativistic electrons. 

In total there are seven free parameters. These are the 
blast-wave energy S52 in units of 10®^ erg, the number den- 
sity no at 10^^ cm (regardless of the density structure), the 
index p of the electron power-law distribution, the index k of 
the density distribution of the matter surrounding the GRB, 
the fraction ^ of accelerated electrons, and te and cb denot- 
ing the fractions of internal energy carried by the relativistic 
electrons and magnetic field, respectively. In practice, due to 
a degeneracy of this model (Eicliler & Waxman 2005) a value 
for one of these parameters has to be assumed in order to 
uniquely determine the others. In this work we 'break' the 
degeneracy by assuming f = 1 in all runs, unless otherwise 
stated. 

The flux prescriptions are based on analytic calculations 
of flux scalings during the relativistic (Blandford & McKee 
1976) and Newtonian (Sedov 1959; Taylor 1950) phase of 
the blast-wave dynamics. In these two dynamical regimes 
the flux at every power-law segment of the spectrum has 
been calibrated in terms of p and k. Several hydrodynamic 
simulations of the afterglow dynamics were run and subse- 
quently post-processed using a radiative-transfer code (van 
Eerten & Wijers 2009; van Ecrtcn et al. 2010a). The cal- 
ibration was carried out by matching analytic expressions 
for the flux scalings to these numerical results. The sharp- 
ness of spectral breaks connecting different power laws of 
the spectrum is also expressed as a function of p and k. The 
transition from the relativistic to the Newtonian solution is 
nicely described as a temporal power-law break between the 
asymptotic behaviour of the critical parameters of the spec- 
trum, namely maximum flux Fm, self-absorption frequency 
Fa and synchrotron characteristic frequency of the lowest- 
energy electrons i^m- It is worth noting that the characteris- 
tics (break time and sharpness) of those temporal breaks are, 
in general, unique for every parameter of the spectrum. This 
emphasizes the advantages of simulation-based flux prescrip- 
tions compared to simple analytic models for the transrela- 
tivistic behaviour of observed afterglows. 
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3.2 The cooling break 

A feature of the synchrotron spectrum not covered in the 
treatment of Leventis et al. (2012) is the cooling break, man- 
ifested as a fourth spectral parameter Vc- Its presence in 
the spectrum, however, might be important, especially for 
observations at optical wavelengths and X-ray energies. For 
that reason all the performed fits have been checked for con- 
sistency by calculating the value of Vc according to formu- 
las available in the literature (e.g. Granot & Sari 2002; van 
Ecrtcn & Wijers 2009) and comparing it to the frequencies 
of the observations. The results of the two aforementioned 
studies are compatible. We have chosen to use those of van 
Eerten & Wijers (2009) due to the fact that a general value 
for k is allowed in their prescriptions. The consistency checks 
have been performed throughout the range of observer times 
covered by the data. A value of Vc greater than the observ- 
ing frequencies implies that cooling has not affected the fits 
and the obtained values for the physical parameters are con- 
sistent with the underlying physical model. To the best of 
our knowledge simulation-based analytic prescriptions for 
beyond the relativistic phase do not exist in the literature. 
That being the case, we have used formulas applicable in 
this phase throughout. This extrapolation provides a lower 
limit on the actual value of because its temporal slope 
in the Newtonian phase is shallower than in the relativis- 
tic (van Eerten et al. 2010a), which is sufficient when Vc is 
found not to interfere with the observing frequencies. 

On the other hand, when the value of is found to be 
lower than - or at about the same levels as - the observ- 
ing frequencies a different approach is necessary in order to 
ffrmly constrain the infiuence of cooling on the data. Our 
fitting code has been expanded to include a prescription for 
the position of fc as a function of time. We have made use 
of the formulas from van Eerten & Wijers (2009) by cal- 
culating of that paper. Formally this expression should 
only apply in the case of slow cooling < Vc)- However, 
it is easy to verify (see also Granot & Sari 2002) that the 
expression for in the case of fast cooling gives a similar 
result within a factor of about 2. A few modifications in 
the prescriptions are then required in order to account for 
the infiuence of cooling in the broadband spectrum. When 
Va. < < or Vm < Va. < fc the ouly modification is that 
of appending another break in the spectrum at the cooling 
frequency, beyond which the spectrum steepens by a half 
(Sari et al. 1998). The formula we have used is 



fobs 



1 + 



?i(a2— aa) 



-1/h 



1 + 



''(03-04) 



-1/r- 



(1) 



The first line in eq. (1) describes the first break of the spec- 
trum at the lowest characteristic frequency, while each fac- 
tor on the second line stands for an extra break at pro- 
gressively higher frequencies. The parameters vo, i>i, v-z and 
s, h, r represent the values of the three critical frequencies 
and the sharpness of the spectral breaks they correspond 
to, respectively, while ai, 02, 03 and a4 are the slopes of the 
four power laws present in a spectrum with three breaks. 
Finally, A is the normalising factor of the spectrum derived 
through modelling of the peak flux Fm- 



When Vjn, Vc < the ordering of I'm and v^. does 
not play a role and one retrieves spectrum 3 of Granot & 
Sari (2002). In that case we have approximated the self- 
absorption frequency with the values applicable to the no- 
cooling case. Similarly, when v^, < Vc < I'm (spectrum 5 of 
Granot & Sari 2002) we have approximated both fm and Vi^, 
with their values in the absence of cooling, while the peak 
flux is attributed to v^. Formally, when < I'c < f'm the 
self-absorption break is split in two break frequencies with 
an extra power-law segment between them that has a slope 
of 11/8. We have neglected that effect and used only one self- 
absorption frequency that has the value of v^x from Leventis 
et al. (2012). This frequency connects power laws of slope 
2 and 1/3. In reality, we have found that most of the time 
best-fit values of the physical parameters imply that these 
approximations arc not used since > v^, Vm- However 
there are instances when this is not the case and we address 
these in more detail in Section 5. 

A last issue that needs to be dealt with when cooling 
infiuences the fits is the application of the relativistic formu- 
las for i^c throughout the range of observer times. To assess 
the validity of this application one needs to estimate the 
duration of the relativistic phase of the afterglow in the ob- 
server frame. In the absence of a detailed description for the 
transrelativistic behaviour of the cooling frequency, the most 
general way to do that is by calculating the observer time 
which corresponds to the transition between the relativistic 
and Newtonian asymptotes, tNR (e.g. Piran 2004; Leventis 
et al. 2012). This calculation has been performed for all sets 
of best-fit parameters and is presented along with our main 
results in Section 4. 



3.3 Fitting procedure 

The fitting method we have used is a x^-minimization algo- 
rithm following the downhill-simplex method combined with 
simulated annealing, as explained in van Eerten et al. (2012). 
The errors for the best-fit parameters of GRB 970508 have 
been determined via a Monte Carlo process. In this anal- 
ysis the values of all data points are perturbed randomly, 
based on their error bars, and a new best-fit set of parame- 
ters is calculated for the synthetic data. For every physical 
scenario (class and constraint) this has been repeated 1000 
times from which 683 best fits were drawn to determine the 
range of the parameters' values at a 68.3%, i.e. Icr, confi- 
dence level. 

The fitted parameters were allowed to vary within the 
following ranges (no in cgs units): 10"'' < E52 < 10"*, 10"^ < 
no < 10^ 2.0 < p < 3.5, lO"'' < cb < 1.0, 10"^ < ee < 1.0, 
-0.5 <k< 2.5. 



4 RESULTS 

For all afterglow data sets, we present three classes of mod- 
els. Each class corresponds to a different assumption (or 
the lack thereof) for the value of k. We have run fits for 
fc = and 2, corresponding to constant density CBM (la- 
belled ISM) and a constant-stellar-wind profile (labelled 
Wind), respectively, and fits where fe is a free parameter. 
For each class, a range of microphysics settings has been 
tested. Namely, we have either allowed for both ee and ee 
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Figure 1. Afterglow of GRB 970508. Best-fit light curves for ISM (solid grey line) and Wind (dashed black line) classes. When 
fe is a free parameter, the Wind scenario is retrieved with high precision. The three radio bands are on top and the rest follow 
in order of increasing frequency, spanning near-infrared, optical, ultraviolet and X-ray energies. Data points before 1.5 days have 
been excluded from the fits but appear in grey in the figure. In all bands, data points have \a errors. Triangles depict upper limits 
at the 2(T level. 



to be free parameters, or connected them through a closure 
relation that effectively reduces them to one free parameter. 
Two options for the closure relation have been explored. On 
the one hand we have imposed equipartition (ee = eb) and 
on the other the 'Medvedev' relation (to = ee; Modvodev 
2006). All other parameters have been kept free at all runs, 
apart from ^ which, for every run, has taken the value of 1. 



4.1 GRB 970508 

We have performed several fits both to the full data set 
and to different subsets (radio only, radio and optical only, 
radio, optical and X rays) of the afterglow observations of 
GRB 970508. Radio data alone do not provide enough in- 
formation to determine simultaneously all the parameters. 
However, when fc is frozen (either in the ISM or the Wind 
scenario) and a microphysics constraint is used, the results 
from fitting the radio only, are fairly similar to those from 
fits to the full data set; all best-fit values of parameters are 
less than 50% off in the Wind class and less than a factor of 
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Table 1. Best-fit parameters, with Icr errors, for GRB 970508 in all classes of models and for all microphysics settings. The fits have 
been performed to data including radio, near-infrared, optical, ultraviolet and X rays, but excluding observations made prior to 1.5 
days after the gamma-ray trigger (see Section 2). The uncertainties in the last row have been calculated for 5 times higher error 
bars. The third column contains the blast- wave energy in units of lO'*'^ erg. The fourth column represents no, the number density at 
a radius of IQ-''^ cm. When fe = 2, no and A, (Chevalier & Li 2000) are related by the formula no ~ 30 A, . For example, the best-fit 
model (equipartition constraint) of the Wind class has A, = 0.243 gcm~^. The last column presents the value of divided by the 
degrees of freedom (dof). 
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2 off in the ISM class. Including X-ray data has almost no 
influence on the inferred values of the physical parameters, 
as the fits arc governed by the combiiiatioii of radio and op- 
tical observations. Nevertheless, -we present results and light 
curves from fits to all bands for completeness. 

In Fig. 1 we present light curves of best-fit models ap- 
plied to the full data set. We have found that the spheri- 
cal model can produce an adequate fit to the data, when 
k = 2. Results for the Wind scenario are almost identical to 
those from fits where fe is a free parameter. Models of the 
ISM class consistently overpredict late radio fiux at 4.86 and 
8.46 GHz. On the other hand, Wind models provide a good 
description at all observer times. In the optical and near- 
infrared bands, the ISM and Wind cases are practically in- 
distinguishable. One common feature of both is the system- 
atic, albeit minor, underprediction of early (< 10 days) flux, 
especially in the R and V bands. This is less pronounced 
in the surrounding K, I and B bands. It is worth noting 
that the X-ray data cannot be fitted by any combination of 
parameters. Along with the fact that the flux drops sharply 
after the first two data points, this hints towards a separate 
origin of the early X-ray flux, for example, inverse Compton 
(e.g. Sari & Esin 2001). Alternatively, the high X-ray flux 
at early times could be due to flaring activity, which is not 
temporally resolved due to the poor coverage. 

In the Wind scenario, all critical frequencies lie below 
the near-infrared. On the other hand, both and z/m pass 
through the radio bands. This is in rough agreement with 
the flndings of Chevalier & Li (2000) and Panaitescu & Ku- 
mar (2002), although we do not conflrm the expectations 
of the former group regarding the passage of Vc from the 
optical. Instead we find that Vc stays below 10^** Hz during 
the observations. In the ISM case we find that starts 
off between the optical and radio and crosses z/a (5 • 10^ Hz) 



at ~ 50 days. We also find that Vc remains between opti- 
cal and X-ray energies throughout, contrary to the results 
of Galama et al. (1998c) and Wijers & Galama (1999) who 
find that Uc crosses the optical frequencies early on. Calcu- 
lation of tNR yields 145 and 180 days, in the best-fit models 
of the ISM and Wind class, respectively. 

In Table 1 we present best-fit parameters, with la er- 
rors, of runs to the full data set. A readily apparent feature 
is the value of k when it is a free parameter, which converges 
to the Wind scenario. Actually, all best-fit values as well as 
the of these two classes are almost identical, regardless 
of the chosen microphysics. From the ISM class only the 
run with no constraints on the microphysics comes close in 
terms of x^, but that model requires a low value for ee and 
high value for ee to work. The energy inferred in this case is 
an order of magnitude higher than the values corresponding 
to the Wind scenario. 

For all classes of models, the best-fit values of x^/dof 
are much higher than 1. This is mainly caused by the notable 
scatter that data in radio, near-infrared and optical bands 
show. The scatter (discussed in Section 2.1) is not refiected 
in the size of the error bars. This is clearly demonstrated 
in the very small uncertainties that the inferred parame- 
ters have, when a microphysics constraint is used. To obtain 
a better measure for the uncertainties when scatter is ac- 
counted for, we have artificially increased the error bars of 
all the data by a factor of 5 and re-calculated them for the 
best-fit model of the k free class. The results are presented 
in the bottom row of Table 1. The choice of the factor is 
motivated by the value of x^/dof ~ 1 that it results in, 
producing a statistically 'good' fit. However, scatter is not 
the only reason for the high values of x^/dof, there are also 
systematic deviations from the data (for example in the R 
and V band during the first 20 days). Therefore, strictly 
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Figure 2. Afterglow of GRB 980703. Best-fit liglit curves for ISM (solid grey line), Wind (dotted black line) and k free (dashed 
black line) classes. Radio bands are shown in the top panel. The lower panel contains neax-infrared, optical and X-ray bands. All 
data were taJjen into account for the light curves we present. In all bands, data points have Icr errors. Triangles depict upper limits 
at the 2(T level. 



speaking, the method of artificially increasing the error bars 
should not be applied to the whole data set. Nevertheless, 
its application results in uncertainties that represent better 
the parameter range allowed by the data and is not used to 
draw any conclusions on the quality of the fits. 

X-ray data show a preference for the ISM class, but 
hardly influence the fit at all, due to the small number of 
data points. We have investigated the dependence of our re- 
sults (especially those for k) on the relative importance of 
X-ray data by increasing the error bars by a factor of 5 in 
all bands, apart from X rays, and recomputing the uncer- 
tainties in the values of the inferred parameters. The best-fit 



results are essentially identical to those of Table 1 for the 
k-free class. The la uncertainties, while larger than those 
presented in the lowest row of Table 1, exclude the ISM 
scenario. 

A discussion of the spectra, dynamics and inferred pa- 
rameters in the Wind scenario (that produces the best fits) 
is presented in Section 5. 

4.2 GRB 980703 

Another well-sampled afterglow that has been extensively 
modelled in the literature is that of GRB 980703. We have 
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Table 2. Best-fit parameters of each class for GRB 980703. For column description see Table 1. 
Class Constraint £52 no p ee Ce X^/dof 



ISM - 17.82 760.3 2.538 10"'^ 1.0 15.56 



Wind - 1.771 14.220 3.865 0.0377 0.133 2 9.090 



kfree Equipartition 2.546 4.265 3.933 0.115 0.115 1.154 8.691 



performed fits to the full data set, from radio to X rays, and 
we have found that no set of parameters can fit the data. 
The Wind model does better than the ISM, but the best fit 
is obtained for fe « 1.15. 



In Fig. 2 we present light curves from best-fit models of 
all classes. In the radio, the ISM model underperforms com- 
pared to the other classes. In the optical and near-infrared 
none of the models scorns to be able to reproduce the data 
adequately, especially in the low-energy bands. X-ray data, 
on the other hand, can only be described within the ISM 
class. From this general picture we can conclude that the 
physical scenario of synchrotron radiation from a spherical 
blast wave is not realistic for this source. 



For every class of models, we have selected those with 
the microphysics settings that produced the best fits and 
present them in Table 2. For the ISM and Wind class, the 
model that performs better is the one with no constraint 
on the microphysics, whereas when k is free, equipartition 
produces the best x^/dof. Fitting the afterglow of this burst 
we have allowed for p to range between 2.0 and 4.0 because 
requiring p < 3.5 results in values on the edge of the pa- 
rameter space. Both the best-fit model of the k free class 
and the one from the Wind class have very high values for p 
(> 3.8). Their x^/dof values are notably better than those 
of the ISM class. The values of t^K are 100, 1310 and 880 
days for the ISM, Wind and k free class, respectively. Due 
to the overall-bad fits to the light curves and the extreme 
best-fit values of p, we consider the values we obtain unre- 
liable. For that reason we have not calculated any errors on 
the derived parameters for this burst. 



It is worth noting the consensus over the outflow ge- 
ometry of GRB 980703. Several studies infer small open- 
ing angles and jet breaks in the timescale of days- weeks 
(Panaitescu & Kumar 2001; Yost et al. 2003; Frail et al. 
2003). In the spherical model, the very fast decays observed 
in the H, J and R bands, can only be explained by very large 
values of p, that result in steep light-curve profiles. However, 
a more natural explanation of the observed slopes would be 
that the edge of the jet has become visible (Rhoads 1999; 
Panaitescu 2005). We therefore regard the results presented 
in this paper implicit confirmation of the jet geometry in 
the outflow of GRB 980703. 



4.3 GRB 070125 

The afterglow of the exceptionally luminous GRB 070125 

was observed in several bands, lasting more than ten days 
in X rays and about a year in the radio. We find that Wind- 
like models provide the best description of the data, but 
with noticeable outliers and with inferred parameters that 
are fairly extreme {E > 10^^ erg, €e = 1, ee < 10~^). We 
consider the results indicative, but by no means conclusive, 
as additional physics (e.g. jets) may be needed to explain the 
deviations and special conditions are required to account for 
the physical parameters we obtain. 

In Fig. 3 we present light curves of the best-fit models 
from each class. Results for the Wind and k free classes are 

similar to each other and differ siginficautly from the ISM 
class in radio, millimetre and X-ray bands, where the former 
perform better. However, late-time behaviour of the data at 
4.86, 8, 46 and 22.5 GHz, as well as millimetre observations 
are hard to explain within any model. In the near-infrared 
and optical bands all classes produce good fits. In the ultra- 
violet, there is a slight underestimation of the flux levels. In 
X rays, only Wind and k free models are able to describe 
the data. 

In agreement with De Cia et al. (2011) we find that Vc 
lies between optical and X-ray energies throughout the ob- 
servations. Chandra et al. (2008), on the other hand, find 
that they can best explain the data when Vc lies below the 
optical. Calculation of fNR yields 80 days in the ISM case 
and ~ 140 days in the other classes. This ensures that the 
relativistic formula for Vc is valid during near-infrared, op- 
tical, ultraviolet and X-ray observations, that last up to 10 
days after the gamma-ray trigger. In all classes, Vm starts 
off bellow the optical and overtakes i-'a at 30 — 80 days. The 
different temporal evolution of i^a makes for the deviations 
in late radio light curves between the ISM class and the 
others. 

In Table 3 we present the values of the inferred parame- 
ters for the best-fit models of each class. Deviations between 
different classes are moderate. Best x'^/dof values are found 
when no assumption for the microphysics is made. This is 
because to explain the data, all models require a high value 
for ee and a very low one for es- Values of the parameters 
when k is free (model with the best x^/dof) are closer to 
those from the Wind scenario, without, however, matching 
them. The inferred energies arc high in all cases, as are the 
values for p. Given the imperfect fits and the extreme pa- 
rameters we infer, we have not calculated errors for their 
values. 
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Figure 3. Afterglow of GRB 070125. Best-fit light curves for ISM (solid grey line), Wind (dotted black line) and k 
free (dashed bla<;k line) classes. Radio and millimetre bands are on top. The lower panel shows near-infraxed, optical, 
ultraviolet and X-ray bands. Data points before 1.5 days have been excluded from the fits but appear in grey colour 
in the figure. In all bands, data points have 1<t errors. 
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Table 3. Best-fit parameters of each class for GRB 070125. For column description see Table 1. 



Class 


Constraint 


E52 


no 


P 






ee 


k 




ISM 




6.968 


1.053-103 


3.203 


1.12-10- 


-6 


1.0 





14.30 


Wind 




11.85 


1.478-103 


2.717 


1.59-10" 


-6 


1.0 


2 


10.80 


k free 




15.32 


3.062-103 


2.831 


5.19-10- 


-7 


1.0 


1.670 


10.49 



5 THE CURIOUS CASE OF GRB 970508 

GRB 970508 is a unique burst in many ways. Its afterglow 
was only the second ever observed and despite the multi- 
frequency monitoring, in some bands over the period of sev- 
eral months, the inferred physical parameters vary widely 
between difeerent authors (Wijers & Galama 1999; Chevalier 

6 Li 2000; Frail ot al. 2000; Panaitoscu & Kumar 2002; Yost 
et al. 2003; Berger et al. 2004). From our sample, the fits to 
GRB 970508 are deemed the most reliable and the most suc- 
cessful, despite the higher values for x^/dof. There are two 
basic reasons for this. The first is the overall behaviour of 
model light curves that successfully reproduce the trends of 
the data at all well-probed wavelengths. The second reason 
is the stability and convergence that the fits to GRB 970508 
show, especially in the Wind scenario, but also when fc is a 
free parameter and a constraint on the microphysics is im- 
posed. If no constraint is placed on the microphysics and 
A: is a free parameter, we cannot discern between the ISM 
and the Wind scenario. However, once either equipartition 
or the Medvedev formula are used, the results clearly favour 
a wind-type CBM. In this Section we present an analysis of 
the physics implied by the best-fit parameters we obtain for 
GRB 970508 and compare those to values inferred by other 
authors. 



5.1 Microphysics 

At first glance, the best-fit values presented in Table 1 re- 
veal an issue concerning the microphysics of the blast-wave, 
namely, the sum of cb and ee is greater than 1. In fact, in 
order for the outflow to be adiabatic, as assumed by the 
model, at least one of these parameters has to be much 
smaller than 1. A low value of ee ensures that most of the 
energy remains in the blast wave, even if the electrons ra- 
diate efficiently, while a low value of eb moderates the en- 
ergy losses of the electron population. The degeneracy of 
the theoretical model (Eichler & Waxman 2005) which has 
prompted us to freeze ^ = 1 during the fit process, can be 
used to solve this issue. The net effect of this degeneracy is 
that a set of parameters — f^^ E52, n'o — f^^ no, e'e = 
/£e, tfl — /^B, = / C produce the same spectrum as the 
unprimed ones, regardless of the value of (the positive num- 
ber) /. Therefore, the inconsistency implied by the high val- 
ues of eb and to may be seen as evidence that < 1, which 
means that not all electrons are accelerated at the shock. 
Consequently, the values for E52 and no presented in Table 
1 should be viewed as lower limits, whereas those for cb and 
Ce as upper limits. 

Another notable feature of the results for the micro- 



physics in the Wind scenario is that we can not conclusively 
distinguish between the three possibilities (no constraint, 
equipartition, Medvedev relation). Equipartition settings 
seem to be marginally favoured by the better x^/dof val- 
ues these models have, but the Medvedev relation cannot 
be ruled out. The ambiguity of our results is mainly caused 
by the relatively high values that both es and te have. We 
have run fits where ^ was frozen at 0.1 and 0.01 and moni- 
tored the behaviour of the two former quantities. They were 
found to be approximately equal to each other and always 
(as did i?52 and uq) followed the scalings implied by the de- 
generacy relations. This confirms energy equipartition be- 
tween power-law electrons and magnetic field, which is also 
suggested by x^/dof values. 

5.2 Spectra 

In the Wind scenario the synchrotron spectrum starts off 
at 1.5 days exhibiting fast cooling (Sari et al. 1998) with 
the critical frequencies having the following values: Va. = 
1.2 - 10^° Hz, i/c = 1.1 • 10^^ Hz, i/m = 1-4 - 10^3^2. At 5 
days, overtakes Vc, causing the wiggle in the radio light 
curves of the model (see Fig. 1). The flux at the highest- 
frequency power law of the spectrum (where near-infrared, 
optical, ultraviolet and X-ray data lie) is independent of the 
ordering of critical frequencies and therefore no feature is 
observed in those bands during the spectral transition. Af- 
ter 5 days the spectrum settles into the slow-cooling regime. 
During the fast-cooling phase (i.e. before 5 days), almost all 
available data lie above Vm and Vc, there are hardly any sig- 
nificant radio observations during that time. Therefore, our 
approximations for when i^a < z^c < I'm have a negligible 
effect on the fits. Moreover, given that the values of Vm and 
i/c are largely independent of their ordering in the spectrum 
(Granot & Sari 2002; van Eerten & Wijers 2009), the valid- 
ity of our approach towards optical data is ensured. From 5 
days onwards, no approximation is made for the value of any 
of the critical frequencies and the model assumes its most 
accurate form. 

It is worth noting that the best-fit spectra naturally 
explain the spectral evolution (at ~ 100 days) depicted in 
Fig. 5 of Frail et al. (2000), due to the passage of u^- In 
the ISM case crosses the radio earlier, at around 45 
days, something excluded by the data. Frail et al. (2000) 
also find — i GHz at seven days, whereas in our best fit 
= 5.5 GHz, at the same observer time. We consider the 
difference negligible, especially considering the strong varia^- 
tion the light curves show around those observer times, due 
to scintillation. This can be verified by inspection of Fig. 4 
of Frail et al. (2000), where the spectral index between 4.86 



© 2012 RAS, MNRAS 000, 1-15 



Fits of a spherical model to GRB afterglows 11 



and 8.46 GHz varies between 0.4 and 1.6 within the first 
two weeks. We have also checked the claim of Galama et al. 
(1998c) who suggest that Vc is observed to pass through the 
near-infrared bands at ~ 10 days. When all the available 
data from several different bands (if, 7, J?, V) are taken into 
account, we find that the spectral index starts off (at ~ 2 
days) having values consistent with late time observations, 
thus showing no evidence of spectral evolution. 

In the best-fit model of the Wind scenario, lies be- 
low the optical bands throughout the duration of optical 
observations. Therefore, its exact value is important at all 
observer times. As mentioned in Section 4.1, calculation of 
iNR yields ~ 180 days. This implies that the values of 
during late near-infrared and optical observations (extend- 
ing up to ~ 200 days in the R band) should be mildly af- 
fected by the transition towards the Newtonian dynamical 
phase. Since Uc is not included in the treatment of Leventis 
ct al. (2012) we do not have a description of the transition 
for this critical frequency, at least not at the level of accu- 
racy that we do for the others. Assuming that the trans- 
rclativistic behaviour of Vc is similar to those of the other 
spectral parameters (smoothly broken power-law) and that 
the break is centered around tNR, we have explored vari- 
ous sharpnesses for that transition and found that the fit 
results remain consistent. The only parameter that changes 
noticeably is p which grows from 2.28 to about 2.34 when 
the transrelativistic evolution of Uc is taken into account. 
Having established that this evolution does not affect the 
inferred parameters, the results we present in Table 1 are 
obtained using the relativistic formula for only. 

5.3 Transrelativistic phase 

First noticed in van Ecrton ct al. (2010a) and subsequently 
quantified in Leventis et al. (2012), the duration of the trans- 
relativistic phase of a spherical outflow in the observer frame 
can be long (this also holds in the case of a jetted outflow; 
Zhang & MacFadycn 2009). The near- infrared and optical 
light curves in Fig. 1 show strong deviations from the ultra- 
relativistic behaviour already at a few tens of days, in ob- 
server time. Their progressive steepening is caused entirely 
by the dynamics slowly adjusting to the Sedov- Taylor solu- 
tion, as there is no critical frequency crossing these bands. 
The effect is similar in the radio, but less pronounced due 
to the simultaneous spectral evolution. 

Deviations of the observed radio light curves from the 
relativistic scalings at timescales of several weeks prompted 
Waxman et al. (1998) to propose a jetted outflow for 
GRB 970508. In this paper we demonstrate how accurate 
modelling of the transrelativistic phase can account for the 
deviations from the ultrarelativistic scalings at observer 
times <IC tNR. This implies that a similar trend may hold 
for at least some other GRB afterglows, the temporal evo- 
lution of which has been interpreted as evidence for a jet 
break. 

Differentiating between jet breaks and the transition to 
the non-rclativistic phase is important, as it directly affects 
the inferred geometry and energetics of GRB outflows. There 
are two main quantities that can serve as diagnostics for this 
differentiation. The first is the change of the temporal index 
of the flux. In the case of a jet break, a decrease in the 
value of the temporal index is expected, mainly due to the 



missing-fiux effect (Panaitescu et al. 1998) that arises when 
the edges of the jet become visible. On the other hand, the 
change in the temporal slope, as the outflow approaches the 
transition to non-relativistic velocities, may be positive or 
negative, is a function of k and (depending on the spectral 
regime) p, and is known from theory (e.g. van Eerten et al. 
2010a). A second diagnostic is the duration (smoothness) of 
the change in the temporal index. In the case of a jet break 
the transition lasts from factors of few (ISM environment) 
up to a decade (wind environment) in observer time (De 
CoUe et al. 2012b; van Eerten & MacFadyen 2012a), whereas 
the typical duration of the transrelativistic regime in the 
case of spherical outflows is a few decades (Leventis et al. 
2012). The picture is slightly more complicated in the case 
of a jet break observed off axis. In that case the jet-break 
transition is effectively stretched and postponed (in reality 
it splits in two). Good coverage is then critical to discern 
between the different interpretations. 

5.4 Comparison to previous work 

Several broadband flts to the afterglow of GRB 970508 have 
been performed and presented in the literature (Wijers 
& Galama 1999; Chevalier & Li 2000; Yost et al. 2003; 
Panaitescu & Kumar 2002). Others have flt only late-time 
radio data (Frail et al. 2000; Bcrger et al. 2004), while Star- 
ling et al. (2008) have flt only the slopes of light curves and 
spectra to infer values for p and k. Most of these studies 
assume or flnd that an ISM scenario flts the data better, 
apart from Chevalier & Li (2000) and Panaitescu & Ku- 
mar (2002) who favour the Wind case. In this study we 
have presented a detailed investigation of both density struc- 
tures that clearly favours a stellar- wind CBM. In addition we 
demonstrate how models with no assumption on the slope 
of the CBM converge to the Wind scenario. Interestingly, 
Starling et al. (2008) find that, in their sample, four out of 
five afterglows with well-constrained values for k suggest the 
same. In that study the density structure of GRB 970508 is 
poorly constrained. 

There seems to be more agreement on the geometry 
of the outflow of GRB 970508. Most studies (also Rhoads 
1999) do not need to invoke a jet, while those that do in- 
fer a jet geometry, usually find large half-opening angles: 
18° (Panaitescu & Kumar 2002), 30° (Frafl ct al. 2000), 50° 
(Yost et al. 2003). We find that the spherical model pro- 
vides a good description of the data, capturing the trends 
of the light curves at different wavelengths for more than 
two orders of magnitude in observer time. We argue that 
GRB 970508 may have indeed originated from an almost 
spherical outflow. The energy of the prompt emission is es- 
timated around 5 • 10*^ erg, if isotropic (Bloom et al. 2001). 
Although on the high side, this value is not unreasonable 
(e.g. Metzgcr et al. 2011). 

In terms of the whole set of fitted parameters, our 
results are similar to those of Chevalier & Li (2000) and 
Panaitescu & Kumar (2002). Given the uncertainties in the 
last row of Table 1, their best-fit values are within, or just 
outside the allowed range of our results. We find moderately 
higher values for £e and es than both studies, but these val- 
ues are effectively upper limits. Lowering ^ to 0.3 results in 
E52 ^ 0.4, At ~ 0.73g cm"\ ee ~ ee 0.19, while the 
value for p remains the same, 2.28. None of these parame- 
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ters are more than a factor of three off compared to both 
aforementioned studies (note, however, the inference of a 
jot from Panaitcscu &: Kumar 2002). It is worth mention- 
ing that the blast-wave energy inferred through modelling 
of the afterglow radiation is very similar to the radiative 
output of the prompt emission. This result holds regardless 
of the outflow geometry and implies a very high cfBcicncy 
of the gamma-ray radiation from the main burst. However, 
given the fast cooling at early times, adiabatic evolution of 
the blast wave demands ee 1. For ^ < 0.17, both ee and ee 
are smaller than 10%. The corresponding blast wave energy 
becomes > 8 • 10^^ erg, which reduces the efficiency of the 
prompt emission below 40%. 

6 DISCUSSION 

In this Section wc discuss the implications of our results for 
the properties of GRB outflows and afterglow fitting. 

6.1 Collimation of GRB outflows 

We have demonstrated how a spherical outflow can account 

for the observations of GRB 970508 and how it fails in the 
case of GRB 980703. The former afterglow has often been 
successfully modelled both with a spherical and a coUimatcd 
outflow. For GRB 980703 a jet is invariably inferred and in 
this research, similarly to Frail et al. (2003), we find that a 
spherical model cannot provide an adequate description to 
the data under any combination of physical parameters. 

The degree of collimation in the case of GRB 070125 
is less clear. On the one hand, there are only a few studies 
of the afterglow radiation and only one of them performs 
broadband (radio to X rays) fitting (Chandra et al. 2008). 
On the other hand, our results provide a satisfactory de- 
scription of most of the broadband data, apart from late 
time behaviour in the radio, when an additional component 
is observed in the light curves. This component, however, 
cannot be explained in the jotted model of Chandra et al. 
(2008) either. Moreover, they propose that the X-ray fiux 
is dominated by inverse Compton, in order to explain what 
seems to be a chromatic break in the optical and X-ray light 
curves ( at about 4 and 10 days, respectively). However, the 
claim for a jet-break in the optical is based only on two data 
points (one in the / and one in the R band), while in the X 
rays it is only based on one (see Fig. 3) . We find that a spher- 
ical model offers a similar level of accuracy, without the need 
to invoke a jet or other radiation mechanisms beyond syn- 
chrotron. However, the parameters we obtain are extreme, 
both on the microphysics side, but also in the total energy 
budget they imply (> 10^* erg). We therefore consider it 
likely that the model wc have used lacks some physics, which 
at least in some bands and observer times 'drives' the radi- 
ated spectrum. That extra physics could be a jetted outfiow, 
but the evidence from previous studies combined with our 
findings is not conclusive. 

In this study wc cannot quantify the opening angle of 
jets, in cases that one is inferred. We can, however, qualify 
afterglows as spherical by successfully fitting their broad- 
band data set. This has been the case for GRB 970508 and 
we consider this a clear demonstration of the diversity in 
the geometry of GRB outflows. This is in accordance with 



searches for jet breaks in large samples of afterglow ob- 
servations that fail to clearly identify a jet break in more 
than half of the sources (Kocevski & Butler 2008; Racusin 
et al. 2009). However, in the coUapsar model (MacFadyen 
& Woosley 1999) for long GRBs, it is Ukely that outflows 
are still coUimated right after breaking out of the stellar 
envelope (e.g. Morsony et al. 2007). If the opening angle 
of the jet is large, the light curves will exhibit deviations 
from spherical symmetry during the transrelativistic phase. 
In such a quasi-spherical scenario, the observational signa^ 
tures of decoUimation might be weak and the expected dif- 
ferences from a perfectly spherical outfiow have, to the best 
of our knowledge, not been explored in the literature. In the 
case of GRB 970508, ^nr = 180 days in the Wind class, and 
the decoUimation should occur on similar timescales, close 
to the end of data sampling. Therefore, if the outflow of 
GRB 970508 had a very large opening angle, our fits will be 
dominated by observations of an almost conical flow. The in- 
ferred energy would then be the isotropic equivalent of the 
real energy content of the blast wave, which is lower only 
by a factor of order unity. The best-fit values of p and k 
are inferred by the slopes of spectra and light curves which, 
at least for the best part of the observations, are not influ- 
enced by effects caused by a possible quasi-spherical geome- 
try. Therefore, we would not expect our general conclusions 
concerning the slope of the CBM to be significantly affected 
by such a scenario. 

Quantifying the distribution of jet opening angles is 
not an easy task, especially considering the inadequate (for 
broadband modelling) coverage that a large fraction of after- 
glows have. On the observational side, Curran et al. (2008) 
have shown that jet breaks may be misidentified as sin- 
gle power laws, due to data-analysis effects. Moreover, van 
Eerten et al. (2010b) have shown that a moderately off-axis 
viewing angle (but smaller than the jet semi-opening an- 
gle) can 'mask' the appearance of a jet-break. If jets are 
present, observing them off a:xis should happen more often 
than not. Therefore, this is an important effect that should 
be taken into account in the model fits. Another issue that 
needs to be better understood is the early (10^ — 10^ s) af- 
terglow behaviour which in a large fraction of bursts sug- 
gests some form of energy injection, continuous or irregular 
(Nousek et al. 2006; Panaitescu & Vestrand 2011). This may 
affect the overall dynamics of the outflow but also result in 
misinterpreting a potentially coincident jet break (Racusin 
et al. 2009). Thus, connecting the dynamics of the early af- 
terglow with the more regular behaviour observed at larger 
timescales is essential to uncover evidence for jets that may 
not be in the form of the canonical achromatic jet break. 

6.2 The immediate environments of gamma-ray 
bursts 

In this work we have treated the density structure of the 
CBM as a free parameter (fc), assuming that a constant 
power law applies. Out of the three data sets we studied, 
one (GRB 970508) showed convergence to a constant stellar 
wind, represented by A: = 2. The best fit to GRB 980703 is 
obtained for k = 1.154. Lastly, for GRB 070125 the best-fit 
value of A; = 1.67, which is closer to that of a constant stel- 
lar wind than homogeneous CBM. For all data sets. Wind 
environments produce better fits than the ISM class. 
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Similarly, however, to the discussion on the geometry of 
the outflows, GRB 970508 is the only one with reliable re- 
sults. In both GRB 980703 and GRB 070125 large values for 
both p and k are needed to best describe the data within the 
spherical model, the applicability of which is at least doubt- 
ful in these cases. For GRB 970508, the value of A* implies 
that the inferred density profile corresponds to a constant 
mass-loss rate of 2.4 ■ 10^*' i^^"' Afe/yr, for a wind veloc- 
ity of l,000kms~^. Interpreting the adiabatic condition as 
^ < 0.17, we find M > 1.5-10"^ Mq/jt, which implies a rel- 
atively massive Wolf-Rayet star towards the end of its life 
(Chevalier & Li 1999). 

Several studies have fit individual bursts and found or 
assumed a homogeneous density structure for the CBM. 
When the fits are compared against those with stellar-wind 
CBM the results are often ambiguous (e.g. Frail et al. 2003; 
Chandra et al. 2008), while in some cases the Wind scenario 
seems to be favoured (Chevalier & Li 2000; Panaitescu & 
Kumar 2002). On the theoretical side, van Eerten & Mac- 
Fadyen (2012a) have shown that the majority of Swift post 
jet break slopes are not reconcilable with a constant den- 
sity CBM, if late energy injection and viewing angle do not 
significantly affect the observations. Instead, the observed 
slopes suggest a wind-type environment for the CBM. Star- 
ling et al. (2008) have studied a sample of 10 Beppo-SAX 
afterglows and found that the majority of the data sets that 
were sufficient to constrain the value of k implied a stellar- 
wind CBM. However, half of them have error bars that allow 
for a wide range for k. Curran et al. (2009) have performed 
a similar study using Swift bursts and find a division in the 
sample between constant and wind-like profiles. It seems, 
therefore, likely that the density structure of the CBM in 
GRBs is diverse, similar to the geometric characteristics of 
their outflows. However, this docs not necessarily translate 
to diversity of the progenitors as well, because in the col- 
lapsar model (Woosley 1993; MacFadyen & Woosley 1999) 
the CBM of a large fraction of long GRBs is modified by 
multiple stellar winds from the neighbouring stars (Mimica 
& Giannios 2011). 

6.3 Model constraints 

For all the afterglows we studied, we have found that a multi- 
frequency data set is more suitable for fitting all the parame- 
ters at once. This has led to the expansion of the model with 
the inclusion of the cooling frequency of the synchrotron 
spectrum, Pc- However, even when radio to X-ray data are 
fitted and all details of the spectrum are taken into account, 
setting k a free parameter results in large uncertainties, if 
no assumption for the microphysics is made. This is mani- 
fested in the large errors for the best-fit values of physical 
parameters in the case of GRB 970508 (see row 7 of Table 
1)- 

When k is free and no microphysics assumption is made, 
the number of fitted parameters is six, equal to the max- 
imum number of constraints we can have from the light 
curves - four from the positions of the critical frequencies 
and the value of Fm , plus two more from the slopes of spectra 
and light curves. However, our results imply that not all of 
these constraints are efficiently used during the fitting pro- 
cess. This means that the effects of two or more of the con- 
straints cannot be separated, leading to a case-specific de- 



generacy. In the case of GRB 970508, for the best-fit model, 
both Vm and Vc lie between radio and near-infrared bands 
for the best part of the observations (vm stays above the 
radio bands for about 100 days). Therefore, their positions 
are not independently constrained by the data, leading to 
a wide range of possible values when all six parameters are 
simultaneously fitted. 

An interesting feature of the prescriptions we have used 
is the inclusion of ^ as a parameter. Due to the degeneracy of 
the model, the presence of ^ is not necessary per se. One can 
imagine a situation where a range in the allowed values for ^ 
is reflected in the adjustment of the ranges of the other pa- 
rameters. For example, by assuming that ^ = 1 and allowing 
ee and Ce to obtain values > 1 during fitting, one accounts 
for the possibility of ^ being smaller than those two param- 
eters, while all of them are smaller than unity. However, the 
inclusion of ^ in the model demonstrates these situations 
more clearly. In the results we obtain for GRB 970508 it was 
not initially possible to discern between the Medvedev con- 
straint and the equipartition constraint for the microphysics 
due to the high values of both eb and ec, that, within the 
uncertainties, extend to the upper limit of the allowed range. 
By freezing ^ at values much lower than 1, we have excluded 
the presence of better fits in which eb > C and/or Ce > 
and confirmed that energy equipartition between power-law 
electrons and magnetic field describes better the afterglow 
observations of GRB 970508. 



7 CONCLUSIONS 

We have performed broadband fits of three afterglow data 
sets using accurate analytic fiux prescriptions applicable to 
spherical outfiows. We have shown that GRB 970508 is suc- 
cessfully fit by a spherical model. The fits fail in the case of 
GRB 980703 and GRB 070125 at varying degrees, implying 
that these sources may be indeed related to jetted outfiows. 
This is supported by extensive modelling of the former and 
the extremely high isotropic energy inferred for the latter. 

For GRB 970508 we find that the best-fit value for k 
is practically 2, strongly suggesting a stellar- wind environ- 
ment. Fits to GRB 970508 also show strong evidence for a 
population of electrons that is not accelerated at the forward 
shock. The implied values for the microphysics parameters, 
te and es, suggest that they are close to equipartition. 

Modelling of GRB 970508 illustrates how an accurate 
spherical model accounts for the progressive deviations of 
light curves from the ultrarelativistic scalings at tobs ^ ^nr. 
This feature had been previously interpreted as a jet break 
in the context of simpler models, but emerges naturally from 
precise calculations of dynamics and spectra in the spherical 
scenario. Therefore, we consider it possible that similar fea- 
tures in the data sets of other afterglows have been misinter- 
preted as jet breaks, in the absence of detailed calculations 
for the spherical case. 
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